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This article reports a measurement of the production cross section of prompt isolated photon 
pairs in proton-antiproton collisions at -^/s — 1.96 TeV using the CDF II detector at the Fermilab 
Tevatron collider. The data correspond to an integrated Inminosity of 5.36 fb~^. The cross section 
is presented as a fnnction of kinematic variables sensitive to the reaction mechanisms. The resnlts 
are compared with three perturbative QCD calculations: (1) a leading order parton shower Monte 
Carlo, (2) a fixed next-to-leading order calculation and (3) a next-to-leading order/next-to-next- 
to- leading- log resummed calculation. The comparisons show that, within their known limitations, 
all calculations predict the main features of the data, but no calculation adequately describes all 
aspects of the data. 
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The measurement of the production cross section of 
two energetic isolated central photons (diphotons) in high 
energy hadron collisions is important for testing standard 
model (SM) predictions in the domain of searches for 
undiscovered particles and new physics. Understanding 
the reaction mechanisms in the complicated environment 
formed in such collisions is a challenge for perturbative 
Quantum ChromoDynamics (pQCD) calculations. Pho- 
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tons orginating from hard collisions of hadrons ( “direct” 
or “prompt” photons) are an ideal probe for testing these 
calculations because they do not interact with other final 
state particles, and their energies and directions can be 
measured with high precision in modern electromagnetic 
calorimeters. Prompt diphoton production creates an ir- 
reducible background to the diphoton decay channel of 
proposed new particles, such as low mass Higgs bosons 
or Randall- Sundrum gravitons in models of extra spatial 
dimensions 0,0. An improved knowledge of the SM 
background will help the development of more powerful 
search strategies for these particles. 

The basic mechanisms of prompt diphoton produc- 
tion in hadron collisions are quark- antiquark annihilation 
qq — t 77, quark-gluon scattering gq — J. 779, and gluon- 
gluon fusion gg — >■ 77. The respective basic diagrams 
are shown in Fig. [TJ At the Tevatron, the dominant 
mechanism is quark-antiquark annihilation. In quark- 
gluon scattering, most of the time at least one of the 
two photons is emitted almost parallel to the scattered 
quark. Contributions from this mechanism are therefore 
suppressed by requiring isolated prompt photons. Each 
mechanism can be modeled by calculating the respective 
matrix element for the specific event kinematics. Ma- 
trix element calculations of Leading Order (LO) in the 
strong coupling are relatively simple and are thus imple- 
mented in advanced parton shower Monte Carlo (MC) 
event generators [10 , which allow for gluon and photon 
radiation as well as multiple interactions in the collid- 
ing beams. By including radiation before and after the 
hard scattering, parton shower generators take into ac- 
count soft gluon and photon emissions, thus resulting in 
an effective resummation of all of the Leading Logarith- 
mic (LL) terms in the cross section to all orders of the 
strong and electromagnetic couplings constants. Next- 
to-Leading Order (NLO) calculations 0-0 additionally 
include one-loop corrections at the cost of not featuring 
realistic multi-particle event representations as the LO 
generators do. Recent NLO calculations include an ana- 
lytical resummation of the cross section for initial-state 
gluon radiation to all orders in the strong coupling con- 
stant 0, reaching a higher logarithmic accuracy than 
in the parton shower Monte Carlo generators. By this 
method, all soft gluon emissions in the initial state are 
taken into account, and reliable predictions for the low 
diphoton transverse momentum region are possible. A 
fixed-order NLO calculation implemented by the diphox 
program 0 also accounts for the case where a final state 
quark loses almost all of its energy to the photon de- 
tected in the event 0 . This process is called “fragmen- 
tation” and, in contrast to final state photon radiation 
in parton showering, it involves non-perturbative calcula- 
tions. One or both photons in the event may come from 
fragmentation. The case where both photons come from 
fragmentation of a single quark is also possible, but is 
not included in calculations, as in this case the photons 
are nearly collinear and non-isolated most of the time. 

The prompt diphoton cross section has been previously 
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FIG. 1: Basic diagrams for prompt diphoton production: (a- 
b) direct, (c-d) one-photon radiation from an initial- (ISR) 
or final-state quark (FSR), (e) fragmentation where one pho- 
ton is emitted along the direction of a final-state quark tak- 
ing almost all of its energy. The symbol ® denotes the 
non-perturbative mechanism of the fragmentation process 
(FRAG). 



measured by the CDF Collaboration using 200 pb ^ of 
data but the large statistical uncertainties did not 
allow for a precise comparison with theoretical calcula- 
tions. The nearly 30 -times larger CDF II data set cur- 
rently available presents an opportunity to significantly 
extend the kinematic range and perform a detailed study 
of diphoton kinematic distributions. A recent measure- 
ment of the diphoton cross section using 4.2 fb“^ has 
been reported by the DO Collaboration The re- 

ported differential cross sections were only partly repro- 
duced by theoretical calculations 000, althormh the 
discrepancies between the NLO calculations 0 , 0 and 
the data were less important in kinematic regions where 
the Higgs boson or new heavy particles are expected. 

This article is organized as follows. An overview of 
the detector is given in Section |B 1 The event selection is 
presented in Section llBl Section HVl deals with extracting 
the cross section from the selected diphoton sample. The 
results are presented and discussed in Section El The 
conclusions are given in Section lVTl Appendix lAlexplains 
details of the non-prompt photon subtraction technique 
introduced in Section HVl Finally, tables of the measured 
cross section, differential in various kinematic quantities, 
are given in Appendix [B] 



II. DETECTOR OVERVIEW 

The CDF II detector is a cylindrically-symmetric ap- 
paratus dll designed to study pp collisions at the Fermi- 
lab Tevatron. The detector has been described in detail 
elsewhere [l0|; only the detector components that are 
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relevant to this analysis are briefly discussed here. The 
magnetic spectrometer consists of tracking devices inside 
a 3-m diameter, 5-m long superconducting solenoid mag- 
net which provides an axial magnetic field of 1.4 T. A 
set of silicon microstrip detectors (LOO, SVX, and ISL) 
[IMi and a 3.1-m long drift chamber (COT) with 
96 layers of sense wires measure momenta and trajecto- 
ries (tracks) of charged particles in the pseudorapidity 
regions of |? 7|<2 and |? 7|<1 [13, respectively. Surround- 
ing the magnet coil is the projective-tower-geometry sam- 
pling calorimeter, which is used to identify and measure 
the energy and direction of photons, electrons, and jets. 
The calorimeter consists of lead-scintillator electromag- 
netic and iron-scintillator hadron compartments and it 
is divided into a central barrel (|? 7 |< 1 . 1 ) and a pair of 
“end plugs” that cover the region l.l<|?y|<3.6. The cen- 
tral calorimeter is composed of towers with a segmen- 
tation of Ar] X A(j) — 0.1 X 15°. The energy resolution 
of the central electromagnetic calorimeter for electrons 
is a{ET)/ET = 13.5%/v'£^T(GeV) 0 1.5% while 
the energy resolution of the central hadron calorimeter 
for charged pions that do not interact in the electro- 
magnetic section is a{E^)/E-r = 50%/i/ifT(GeV) © 3% 
[l3 • Multiwire proportional chambers with cathode-strip 
readout (the CES system), located at the depth of six 
radiation lengths (near shower maximum) in the central 
electromagnetic calorimeter, are used for identification 
and precise position measurement of photons and elec- 
trons. Gathode strips and anode wires, with a channel 
spacing between 1.5 cm and 2 cm, running along the 
azimuthal (strips) and the beam line (wires) direction 
provide location and two-dimensional profiles of electro- 
magnetic showers. The position resolution of the GES is 
2 mm for a 50 GeV photon. The electromagnetic com- 
partments of the calorimeter are also used to measure the 
arrival time of particles depositing energy in each tower 
M. A system of Cherenkov luminosity counters (CLC) 
2 y, located around the beam pipe and inside the plug 
calorimeters, is used to measure the number of inelastic 
pp collisions per bunch crossing, and thereby the lumi- 
nosity. 

The online event selection at CDE is done by a three- 
level trigger ( 2 ^ system with each level providing a 
rate reduction sufficient to allow for processing at the 
next level with minimal deadtime. Level- 1 uses custom- 
designed hardware to find physics objects based on a 
subset of the detector information. Level-2 does limited 
event reconstruction. Level-3 uses the full detector infor- 
mation and consists of a farm of computers that recon- 
struct the data and apply selection criteria similar to the 
offline requirements. 



III. DATA SELECTION AND EVENT 
RECONSTRUCTION 

Inclusive 77 events are selected online by a three-level 
trigger that requires two isolated electromagnetic (EM) 



clusters with E^>12 GeV (diphoton-12 trigger) or two 
electromagnetic clusters with E^>18 GeV and no isola- 
tion requirement (diphoton-18 trigger). The transverse 
energy of the clusters is calculated with respect to the 
nominal center of the detector at z=0 cm. The trigger 
requirements at each level are briefly described below. 

At Level-1, events having two towers with EM Et>8 
GeV each are required. Eor each trigger tower, the 
amount of energy in the hadronic compartment of the 
calorimeter has to be consistent with that of an 

electromagnetic object. A trigger tower consists of two 
adjacent towers in the same calorimeter wedge, so that 
the granularity is approximately Arj x A(f> ~ 0.2 x 15°. 

The Level-2 requirements are different for the two trig- 
gers. The diphoton-12 trigger selects events if there are 
two isolated seeds with EM Et> 10 GeV each. The isola- 
tion (ISO) energy is calculated as a sum of the transverse 
energy in the towers nearby the seed tower. The ISO en- 
ergy for both photons has to be less than 3 GeV or 15% 
of the seed energy, whatever is larger. The diphoton-18 
trigger requires two towers with EM At >16 GeV each at 
Level- 2. 

Events are fully reconstructed at Level-3. At this level, 
for all photons in both triggers, the energy profile at the 
shower maximum of each photon candidate has to be 
consistent with that of a single photon. The diphoton- 12 
trigger selects events with two isolated photon candidates 
with At >12 GeV. The isolation energy at level-3 is cal- 
culated as the sum of At in all towers (except for photon 
towers) within the cone of AR = y/ {Arj)^ + (A^)^ < 0.4 
centered around the photon candidate. This ISO energy 
has to be less than 2 GeV or 10% of the photon energy, 
whatever is larger. The diphoton- 18 trigger has no iso- 
lation requirement and accepts events with two photon 
candidates with At >18 GeV. Table U gives a summary 
of all trigger requirements for events with EM objects in 
the central calorimeter and with At calculatated with 
respect to the event vertex. 

The triggered 77 candidate events are then subject to 
the offline selection. Each event is required to have two 
central photon candidates inside a well-instrumented re- 
gion of the calorimeter (approximately 0.05<|t7|<1.05) 
with At>17 GeV for one candidate and At>15 GeV 
for the other. This asymmetric cut helps to avoid in- 
stabilities in fixed NLO calculations Q- Photon can- 
didates must satisfy strict (referred to as “tight”) pho- 
ton identification requirements. The EM cluster has to 
be located inside the welEinstrumented region of the 
GES chamber, away from the (()-boundary of a calorime- 
ter tower [ 23 . The energy deposition pattern in both 
transverse profiles at GES has to be consistent with 
that of a single electromagnetic object. The ratio of 
the energy measured in the hadron (HAD) calorime- 
ter to the EM energy, A^^'^/A^^, has to satisfy the 
requirement A^"^^/A'®^<0.055+0.00045x A'*'. To dis- 
tinguish photons from electrons, no high-pT charged- 
particle track should point into the cluster {Ntrack^^ 
with track pt<1-0+0.005x At). The main sources of 
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TABLE I: Summary of the diphoton trigger requirements. 



Trigger Level 


Diphoton- 12 


Diphoton-18 




EM Et > 8 GeV 


same 


Level- 1 


EHAD/eEM ^ Q ^25 


same 




^cluster — 2 


same 




EM Et > 10 GeV 


EM Et > 16 GeV 


Level-2 


EHAD/eEM ^ Q ^25 


same 




< 3 GeV or < 0.15 


not applied 




^cluster — 2 


same 




EM Et > 12 GeV 


EM Et > 18 GeV 


Level-3 


EHAD/eEm ^ Q Q 55 0.00045xE/GeV if E < 200 GeV 


same 




< 2 GeV or E^^°/Et < 0.1 


not applied 




shower profile: Xces < 20 


same 




^cluster — 2 


same 



“fake” photons are energetic 7 r° and rf mesons produced 
in jets. These mesons are usually produced in associa- 
tion with other particles. To reduce this contamination 
from jets, the photon candidate must be isolated in the 
calorimeter. To calculate the calorimeter isolation (cal- 
ISO), the deposited in the calorimeter towers within 
the cone of Ai? < 0.4 around the EM cluster is summed, 
and the Et of the EM cluster is subtracted. Cal-ISO is 
then corrected for the photon’s energy leakage into tow- 
ers in the neighboring wedge and for the contribution 
from multiple interactions in the same bunch crossing 
[ 2 ^. Cal-ISO must be consistent with the amount of en- 
ergy expected from the underlying event (see Table El. 
In addition to the calorimeter isolation, there should be 
no other significant energy (Et of 2"^* CES cluster) de- 
posited in the CES chamber containing the photon can- 
didate. Table HIl provides a summary of the photon iden- 
tification requirements described above. To reduce con- 
tamination due to cosmic-ray, beam-related, and other 
non-collision backgrounds, the event must contain a well- 
reconstructed vertex, formed from tracks, with |z|<60 
cm. If multiple vertices are reconstructed, the vertex 
with the largest^ pt of the associated tracks is selected. 
The transverse energy of the photon candidates is calcu- 
lated with respect to this primary vertex. 

Inclusive 77 events satisfying the above criteria form 
the baseline 77 sample used in the analysis. Due to the 
presence of fakes, this sample consists of real 77 , 
and jet-jet events. (An object misidentified as a pho- 
ton is referred to as a “fake” photon.) Events with one 
or two fake photons are classified as background. The 
baseline signal plus background 77 sample consists of 
roughly 60,000 events in data corresponding to 5.36 fb“^ 
of integrated luminosity. Signal and background samples 
were simulated with the pythia event generator which 
includes simulation of the underlying event and multiple 
hadron interactions, as well as initial- (ISR) and final- 
state radiation (FSR) and a hadronization model of the 
final state partons Q . The pythia events were processed 
through a GEANT-based detector simulation [ 2 ^ and trig- 
ger emulation, followed by the same reconstruction pro- 



gram as that for the data. 



IV. CROSS SECTION MEASUREMENT 

This section describes the steps of the cross section 
measurement. Kinematic variables of interest are his- 
togrammed to measure the corresponding differential 
cross section. The background is subtracted from each 
histogram bin. The signal histograms are normalized to 
the integrated luminosity and to the size of each bin to 
obtain uncorrected differential cross section histograms. 
These are then corrected for the reconstruction efficiency, 
acceptance, and resolution effects. 



A. Background Subtraction 

The fake photon background subtraction is based on 
the use of the track isolation (track-ISO), which is calcu- 
lated as the 'YjPT of tracks with Ai? to the photon <0.4 
and \zyertex~ztrack\<^ Cm. The concept of this technique 
is similar to the one used in the earlier measurement of 
the inclusive photon cross section [l^ . The main idea be- 
hind the method is that true and fake photons have very 
different isolation distributions (see Fig. [2). Therefore, 
one expects different efhciencies for signal (true photons) 
and background (fake photons) for a given isolation cut. 
In a single-photon sample this property can be used to 
extract the number of true photons: 



E e-i — ebjETi) , 

,=i ^s{ETi) — eb{E^i) 

where ti=l if track-ISO<cut and ej=0 if track-ISO>cwf, 
es(ET) is the signal efficiency for track-ISO<cwt, eb[E^) 
is the background efhciency for track-ISO<cut and N is 
the total number of candidate photons in the sample. 
This technique can be generalized in the case of the 77 
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TABLE II: Summary of the standard (“tight”) photon identification requirements for the 77 sample. 



Cuts 


“Tight” photon ID 


Calorimeter fiduciality 


central 


El 


> 15 GeV (U‘ 7 ), > 17 GeV ( 2 “'* 7 ) 


Shower profile in CES: 


<20 


AD j Tp^h/IVl 


<0.055+0.00045xA/GeV 


cal-ISO 


<0.1xAt if At<20 GeV or 
<2.0 GeV-h0.02x(AT - 20 GeV) 


Ntracks in duster 


< 1 


track Pt if Ntracks = 1 


<1.0 GeV+O.OOSxAx 


Et of 2"-“ CES 


^0.14x_£/t if GeV 


cluster 


< 2.4 GeV-tO.OlxAx if At^IS GeV 



sample and is based on a maximum likelihood approach 
which is described in detail in Appendix El 




FIG. 2: The track-ISO distribution in signal (black-solid) and 
background (red-dashed) events. 

As mentioned in Section fllli two types of isolation can 
be defined for central photons: calorimeter and track 
isolation. Cal-ISO is sensitive to the following contribu- 
tions: underlying event (UE), multiple interactions (MI), 
leakage from the photon cluster (LE), and fragmentation 
contribution (FR) from jets (for fakes). Track-ISO, on 
the other hand, is only a measure of UE and FR. There- 
fore, it can potentially offer a better separation between 
true and fake photons. Using track-ISO for the fake 
photon background subtraction also has additional ad- 
vantages that 1ow-Pt tracks are very well measured (un- 
like the calorimeter energy) and jet fragmentation studies 
[ 2 ^ indicate that track observables are well described by 
PYTHIA both for the UE and for the jets. 

To perform the background subtraction, signal (cg) and 
background (ei,) efficiencies are needed for a certain cut 
on track-ISO. The form of Eq. [T] suggests that the best 
accuracy in photon purity can be achieved when the ab- 
solute value of the denominator is maximum. When this 
happens the terms in the sum of Eq. [T] are minimized in 
magnitude and thus the purity is less sensitive to the sta- 



tistical uncertainty of the number of events in the sample. 
Therefore, a scan of Cg — e?, as a function of the track-ISO 
cut is performed using MC samples of true and fake pho- 
tons. The difference peaks at track-ISO~l GeV. The 
threshold of the track-ISO cut for the signal and back- 
ground efficiency functions is thus chosen at 1 GeV. 

The signal track-ISO efficiency is obtained from 
PYTHIA photon-plus-jet samples. The background ef- 
hciency is obtained from pythia dijet samples. All 
PYTHIA samples used in this work are derived from ver- 
sion 6.2.16 of the program using the CTEQ5L PDF 
set for [23 and the “tune A” for UE parameters [^ . 
Background events are filtered out if a detector pho- 
ton is matched to a generator level photon originating 
from quark ISR or FSR. This ensures that the back- 
ground track-ISO efficiency function is obtained for neu- 
tral hadrons (mostly or 77 *^) faking a photon signature. 
Similarly, for the signal events detector photons are re- 
quired to match generator level photons from the hard 
scattering (thus fragmentation photons are removed). 
The signal efficiency Cg and the background efficiency Cb 
are shown in Fig. [3] as functions of the photon Et- Both 
functions are parameterized by a linear combination of 
an exponential and a constant. 

The isolation cones for the two photon candidates are 
not entirely independent. For example, if a particular 
event has a higher (lower) than average underlying event 
activity, then it is likely that both isolation cones will 
simultaneously have more (less) energy. In addition, the 
ordering in Et of the two photons also introduces some 
bias. The signal (for Et<50 GeV) and the background 
track-ISO efficiencies drop with increasing Et (see Fig. 
El). Therefore, the lower Et threshold for the hrst pho- 
ton in the event relative to the second photon implies, on 
average, that the Et of the first photon will be system- 
atically lower than the Et of the second photon, thus in- 
troducing some bias due to the Et dependence of the ef- 
hciencies. This effect is negligible for At>50 GeV, where 
the signal track-ISO efficiency is flat and the background 
is weak, but it becomes significant at low Et- It is a small 
effect in the single photon purity, but it is at least a factor 
of two more important for diphoton events. These cor- 
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:-ISO<l GeV. The shaded area is the total systematic uncer- 



relations must be taken into account when calculating a 
probability of two photon candidates to pass-pass, pass- 
fail, fail-pass or fail-fail the 1 GeV isolation cut described 
above, pythia diphoton events were used to obtain “per 
event” track-ISO efficiencies for these combinations. Cor- 
relations are much less important for events with one or 
two fake photons because they are diluted by a much 
larger contribution from jet fragmentation. 

The systematic uncertainties in the signal and back- 
ground track-ISO efficiencies are estimated and propa- 
gated into the final estimate of photon purity. Correla- 
tions between different sources of systematics are taken 
into account. The following sources of uncertainties are 
considered: 1) mis- modeling of the distribution of the 
number of vertices {Nyx) in MC (pile-up effect); 2) sta- 
tistical uncertainties in the fit parameters; 3) choice of the 
fit function for the efficiency; 4) generator-related data- 
MC differences; 5) effect of the Et threshold for selected 
photon candidates (only for the background). These un- 
certainties are presented in Fig. 2] and discussed below. 

The MC simulation does not describe accurately the 
distribution of the number of reconstructed vertices. This 
effect can be either removed by re-weighting the MC to 
match the data or the associated uncertainty can be as- 
signed for the effect of mis-modeling. The latter approach 
is chosen in this analysis because the track-ISO, to lead- 
ing order, is not sensitive to the presence of multiple in- 
teractions and the effect is very small. This uncertainty 
is conservatively estimated as the difference between the 
extreme cases of track-ISO efhciencies obtained in events 
with Nyx=^ and track-ISO efhciencies obtained in events 
with For the photon energies relevant to this 

analysis, the relative effect is <1% for the signal and 
<3% for the background. 

The ht statistical uncertainties are included in the esti- 
mation of systematic uncertainties. Correlations between 
ht parameters are properly taken into account. The rel- 
ative effect is negligible for the signal and l%-3% for the 
background in the range of photon energies relevant to 
this analysis. 

The default ht function choice is an exponential plus 



a constant term. As one can see from the track-ISO ef- 
hciency plots shown in Fig. [31 the quality of the ht is 
very good. The studies of the efficiency dependence on 
the track-ISO cut indicated that, for some cut values, the 
exponential plus a linear function can be a better ht to 
the signal efficiency. This function was thus chosen as an 
alternative track-ISO efficiency parameterization and the 
difference with the default function was taken as the as- 
sociated uncertainty. The relative effect is <1% for both 
signal and background with At <200 GeV. 

The modeling of both signal and background relies on 
the MC. Therefore, it is necessary to assign a systematic 
uncertainty on possible data-MC differences both for sig- 
nal and background photons. In the case of signal, it is 
necessary to check the modeling of the underlying event 
in the MC. This is done by means of complementary 
cones. The complementary cones are chosen such that 
their axes have the same angle 9 with respect to the beam 
line as the photon candidate and are rotated by ±7 t/ 2 
in 4>. These cones are assumed, on average, to collect 
the same amount of the underlying event as cones of the 
same size around true photons. This assumption is tested 
and confirmed in the MC. It is also checked that comple- 
mentary cones for signal and background look very sim- 
ilar. Finally, the signal track-ISO efficiency is obtained 
from complementary cones in data and signal MC and 
the difference between the two is taken as the associated 
systematic uncertainty in the signal track-ISO efficiency. 
The comparison of track-ISO efficiencies for complemen- 
tary cones in data, MC signal and background is shown 
in Fig. [5] The relative effect is ~3.5% for most of the 
photon energies. This is the largest systematic uncer- 
tainty for the signal track-ISO efficiency. The system- 
atic uncertainty from this source decreases from 5.5% at 
At= 10 GeV to ~3.5% at At= 40 GeV and then it stays 
at roughly the same level for At > 40 GeV. An additional 
uncertainty arises from the fact that the signal efficiency 
is derived from true photons generated only by direct 
gq — >■ 79 , qq — )■ 'yg and gg — >■ yg production, omitting 
photons radiated from initial or final state quarks. Fig. 
|3| shows that the track isolation of photons radiated from 
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plementary cones from data (squares), signal MC (triangles) 
and background MC (circles). 




FIG. 6: The track-ISO distribution in hard scattering and 
ISR (blue) and FSR (black) events. 



final state quarks is somewhat different than that of pho- 
tons produced by hard scattering or radiated from initial 
state quarks. This difference is estimated to have a con- 
stant 2% effect on the signal track-ISO efficiency which 
is added to its total systematic uncertainty. 

The systematic uncertainty due to data-MC differ- 
ences in the track-ISO background efficiency is esti- 
mated by comparing the track-ISO cut efficiency for 
a leading track in dijet events from data and MC. 
This method assumes that jets with a leading neu- 
tral particle (e.g., have the same or very sim- 

ilar fragmentation properties as jets with a leading 
charged particle (e.g., or K=*=). The following pro- 
cedure is applied to both data and MC. Events with 
two well-balanced and back-to-back jets are used, satis- 
fying |ET(jetl)-£;T(jet2)|/[£;T(jetl)-hET(jet2)]<0.3 and 
l^(jetl)— 0(jet2)|>2.7 rad. The event is rejected if there 
is a third jet with Et>0.Ix [ET( jetl)-|-ET(jet2)]. One 
of the jets (a probe jet) is required to be in the cen- 
tral detector region, |? 7 |<I.l, thus matching the pseudo- 
rapidity requirement for photons. In the next step, a 
well-reconstructed track (a probe track) is selected with 
the largest px inside the probe jet, i.e. inside a cone 



of Ai?=0.4 around the jet direction. For this track, an 
analog of the cal-ISO is calculated as the y) Et of all 
towers inside a cone of Ai?=0.4 around the track direc- 
tion. Towers associated with the track (up to 3 towers in 
rj) are excluded from the sum. The cal-ISO for the probe 
track has to satisfy exactly the same requirements as the 
isolation for a photon with Et=Pt- An analog of the 
track-ISO for the probe track is also calculated by fol- 
lowing exactly the same procedure as for photons, with 
the only exception being that the track itself is excluded 
from the sum. Finally, the efficiency of the track-ISO<l 
GeV cut for the probe track is compared in data and 
MC. The observed relative difference of 8%, independent 
of the track px, is taken as an estimate of the systematic 
uncertainty due to data-MC differences in the track-ISO 
background efficiency. 

Finally, the last source of systematic uncertainty in 
the background track-ISO efficiency is associated with 
the choice of an Et threshold for selecting fake photons 
from a particular jet sample. The fake rate for jets is very 
small and, as a consequence, the MC dijet samples do not 
have enough statistics to yield a sufficient number of fake 
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FIG. 7: The estimated signal fraction in the inclusive photon 
data. The shaded area is the total systematic uncertainty in 
the signal fraction. 



photons after the selection cuts. To maximize the statis- 
tics, fake photons are accepted from each dijet sample 
if where px is the parton transverse momentum 

cutoff used in the event generation. Ideally, events with 
i?T>(PT+offset) should have been selected to avoid a bias 
due to the px threshold effect. This is necessary because 
fake photons carry, on the average, only 90% of the en- 
ergy of the original parton. Therefore, the procedure is 
biased toward selecting fakes originating from gluon jets 
produced by radiation, which are not limited by the px 
threshold of hard scattering, and as a consequence to- 
ward lower background efficiencies. To obtain a conser- 
vative estimate of this effect, the threshold was lowered 
even more, thus accepting fake photons with i?x>T'xpx 
where F~0.8-0.9, depending on px- By decreasing the 
threshold, the effect is overestimated, but this gives a 
conservative estimate of the associated uncertainty. The 
total systematic uncertainty of the background track-ISO 
efficiency is at the level of 10%-I2% in the range of pho- 
ton E't from 15 GeV to 200 GeV, the range relevant to 
this analysis. 



The background subtraction procedure has been tested 
with MC signal and fake events as well as with inclusive 
photon data. Tests with MC provide closure checks: the 
returned purity was 100% for signal events and 0% for 
fakes (within the corresponding uncertainties). The esti- 
mated photon purity for inclusive photon data as a func- 
tion of the photon transverse energy is shown in Fig. [7] 
and is similar to the purity obtained in the inclusive pho- 
ton cross section analysis [T^ . The uncertainty in the sig- 
nal fraction of the inclusive photon sample achieved with 
the track-ISO method is between ~11% at low E^ and 
~5% at very high E-[. Fig. [8] shows the estimated pu- 
rity for the diphoton data as a function of the kinematic 
variables defined in Subsection IV Bl 



TABLE III: Summary of the requirements applied to the gen- 
erated MC events to define the denominator of the selection 
efficiency and the meaning of cross section measurement itself, 
is the photon rapidity [l^ . 



Cuts 


Selected 77 events 




> 15 GeV (F‘ 7 ), > 17 GeV (2"^7) 


ll/"l 


< 1.0 for both photons 


Isolation 


<2.0 GeV for both photons 


AR(77) 


>0.4 


No matching requirement between 
generated and reconstructed objects 



B. Event Reconstruction and Selection Efficiency 

The corrections for event reconstruction and selection 
efficiency were derived primarily from pythia diphoton 
MG samples. The numerator of the efficiency is the num- 
ber of events with two photons that pass all of the trigger 
criteria and selection cuts listed in Tables U and n The 
definition of the cross section measurement is determined 
by the definition of the denominator of the efficiency. 
The denominator cuts are summarized in Table Hill This 
work reports a cross section for isolated photons, so the 
selection of denominator events includes isolation. This 
isolation is found by summing over all generated hadrons 
and photons originating from the primary vertex within 
a cone of Ai?=0.4 around each photon. 

For each kinematic quantity, one histogram of the re- 
constructed quantity and one of the quantity derived 
from generator variables are constructed. In the first 
iteration, the efficiency is computed as the ratio of these 
histograms. This ratio also corrects event migration in 
neighboring bins due to finite resolution. Events which 
pass the denominator cuts and have a reconstructed value 
for the histogrammed quantity but not a generator level 
value are assigned the reconstructed values as the best 
approximation to the generator level values. To improve 
accuracy, the efficiency calculation is iterated a second 
time. Once all corrections are applied to the data, includ- 
ing the efficiency, it is the best available representation 
of the true distribution. Then the pythia events are 
reweighted so that the second iteration of the denom- 
inator histogram agrees with the corrected data. The 
purpose is to correct the pythia distribution closer to 
the true distribution, making the efficiency more accu- 
rate. In practice, this does not have a large effect on any 
distribution (see Fig. |9]). A third iteration changes the 
efficiency at the level of 1% or less and, therefore, only 
two iterations are applied. 

The following corrections are applied to the efficiency: 

• — I e+e" events in data and MG are compared 

to derive a correction to the photon ID efficiency 
reported by the MG. The correction is reported as 
a function of and of run periods. The correc- 
tion is weighted by the period luminositites, and 
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FIG. 8: The estimated signal fraction in the diphoton data as 
total systematic uncertainty in the signal fraction. 



the observed Ny^ distributions to find an overall 
multiplicative efficiency correction of 0.967 per pho- 






function of several kinematic variables. The shaded area is the 



ton. In addition, there is some indication of an E^- 
dependence, so the factor 0.967 is allowed to vary 
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linearly up to 1.0 between 40 and 80 GeV, and then 
is held constant at 1.0 above 80 GeV. 

• A small correction is included near the At cut 
threshold due to the trigger turn-on curve. This is 
implemented as a factor of 0.98 at At = 15 GeV, 
going linearly up to 1.0 at At = 18 GeV. Ref. [l^ 
concludes that there is no need for other corrections 
for the trigger. 

• PYTHIA includes the underlying event, but NLO 
calculations do not. This makes the PYTHiA-based 
efficiency correction too large when comparing the 
isolated cross section to NLO predictions. It is too 
large since the UE causes events to be removed from 
the isolated denominator of the efficiency. A cor- 
rection is derived by convoluting the pythia UE 
isolation energy with the diphox energy in the iso- 
lation cone @ . This reduces the probability for the 
DIPHOX event to pass the isolation cuts. This effect 
is measured to be a factor of 0.88 per event which 
is then applied as a correction to the data. 

The efficiency obtained for the kinematic quantities de- 
fined in Subsection IV Bl is shown in Fig. [S] The typical 
efficiency is 40%. 

In addition to the total systematic uncertainty aris- 
ing from the background subtraction, whose details are 
discussed in Subsection IIV A[ the following systematic 
uncertainties are included in the cross section measure- 
ment. 

• The A*^-based efficiency correction has an uncer- 
tainty from several sources, including uncertainty 
in the amount of material leading to conversion 
events, which are rejected, and the difference be- 
tween the electron and photon response to cuts. 
These are summarized as 1.8% below At = 40 
GeV, rising linearly to 3% at 80 GeV and fixed 
above that point. This increase completely cov- 
ers the AT-dependence in the photon ID efficiency 
mentioned above. 

• The photon energy scale is varied and the change 
in the kinematic distribution is reported as an un- 
certainty. For the diphoton mass, the variation is 
0 at At = 40 GeV, rising linearly up to 1.5% at 80 
GeV, then fixed above 80 GeV. These uncertainties 
are based on energy scale studies in the inclusive 
photon cross section measurement [^ . 

• A 3% uncertainty due to trigger efficiency is taken 
from Ref. 0- 

• A 6% uncertainty (3% per photon) for underlying 
event correction is taken from Ref. [1^ . 

• No uncertainty in the acceptance from variations 
in the ISR/FSR model is included since the pri- 
mary mechanism for the effect is extra jets inter- 
fering with isolation. Both the numerator and de- 
nominator photons in the efficiency calculation are 



isolated, therefore the efficiency is immune to this 
effect, to leading order. 

• No uncertainty in the efficiency due to the choice of 
the scale is included because the primary mech- 
anism of this effect is through the boosting of the 
final state. Since the efficiency’s numerator and 
denominator are calculated with full kinematic re- 
quirements, the efficiency is immune to this effect 
to leading order. 

The breakdown of the systematic uncertainties for the 
kinemaric quantities defined in Subsection IV Bl is shown 
in Fig. 1101 In all distributions, the dominant uncertainty 
comes from the background subtraction. The total sys- 
tematic uncertainty is obtained by adding all individual 
components quadratically and averages near 30%. 

C. Corrections and Tests with the e^e~ 

Sample 

The — >■ e+e“ sample is used in this analysis for two 

purposes: 1) to set the energy scale in the GEM calorime- 
ter and 2) to check the overall cross section normaliza- 
tion. The — >■ e+e" data sample is derived from the 

diphoton trigger dataset. The same global event selection 
as for the diphoton sample is applied. Two objects are 
required to pass a “photon-like electron” selection. The 
cuts applied are those of the standard photon selection, 
with modifications to allow for the electron track. The 
modifications are: 

• The number of allowed tracks in the cluster is in- 
creased by 1. 

• The leading track pT cut is applied on the second- 
highest pt track instead of the highest one. 

• The track isolation is corrected by subtracting the 
leading track pT. 

• 0.8 < A/p < 1.2 is required for the energy-to- 
momentum ratio of the leading track. (Events 
which fail this cut also tend to fail the GES 
cut.) 

The electrons are required to match EM objects passing 
the Level- 1 and Level-2 trigger criteria (or trigger sim- 
ulation for the MG). A high-luminosity sample of fully 
simulated and reconstructed pythia — )■ e+e” events 

is used for MG. 

The electromagnetic energy scale is set by tuning the 
reconstructed mass to the world average [s^ in both 
the data and the MG samples. The correction is applied 
as a function of time. It is applied before final event selec- 
tion to account for a few events slightly below the energy 
threshold which the correction pushes above the thresh- 
old. The correction can only have a noticable effect on 
kinematic variables with rapidly falling spectra, related 
to the photon At, such as the diphoton mass. A 1.5% 
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FIG. 9: The estimated efficiency as a function of several kinematic variables. The shaded area is the total systematic uncertainty 
in the efficiency. 



systematic uncertainty due to energy scale is included, as 
mentioned in Subsection IIV Bl 



The ^ e+e- cross section is measured in order to 
check the cross section measurement procedures. This 
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FIG. 10: The estimated systematic uncertainties in the 

measurement tests the trigger efficiency, the ability for 
the MC to predict the event selection efficiency, the effi- 
ciency corrections, and the luminosity. The cross section 






section as a function of several kinematic variables. 



is measured for events with e“'"e“ invariant mass between 
65 and 115 GeV/c^. The “photon-like electron” selection 
is applied and the efficiency from the pythia MC is used. 
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The same photon efficiency corrections as in the dipho- 
ton cross section (see Subsection llVBI) are applied. Since 
the photon-like electron cuts are used, it is assumed that 
the response is similar to photon response. The resulting 
e“'"e cross section is found to be consistent with 
previous dedicated measurements and expectation from 
theory [3lj |. 



V. RESULTS 

This Section presents the results of the cross section 
measurement. A brief description of the theoretical cal- 
culations is given first, then the comparisons for selected 
kinematic variables are shown and discussed. Tables with 
the measured cross section values are given in Appendix 

IS 



A. Theoretical Calculations 

The results of this measurement are compared with 
three theoretical predictions: 

• A calculation using the pythia program Q . This is 
a parton-showering generator which features a re- 
alistic representation of the physics events in terms 
of observable particles. It includes initial and fi- 
nal state radiation and an underlying event model. 
PYTHIA implements a Leading-Order (LO) matrix 
element (ME) for direct diphoton production which 
includes the qq — >■ 77 and gg — >■ 77 LO processes 
described, respectively, by diagrams (a) and (b) of 
Fig. [T] Significant contributions also arise from 
the processes qq 773 (diagrams (c) of Fig. [T]) 
and gq — >■ 77(7 (diagrams (d) of Fig. [1} where 
the second photon is emitted from an initial or fi- 
nal state quark according to the pythia radiation 
model. These contributions were included in the 
calculation by running the program with a filter se- 
lecting diphoton events from inclusive 7-I-X events, 
where X is either a photon or a jet, with an ef- 
ficiency of 0.025%. Fig. [TT] shows the individual 
contributions to the cross section as a function of 
the diphoton invariant mass, transverse momentum 
and azimuthal difference. Initial-state radiation 
(ISR) photons, in particular, produce substantially 
different distributions than ME and final-state ra- 
diation (FSR) photons, having a harder transverse 
momentum spectrum and stronger low~A(/) tail in 
the azimuthal difference spectrum. In leading or- 
der, this can be attributed to the fact that FSR 
occurs in quark-gluon scattering [diagram (d) of 
Fig. [1], whereas ISR occurs both in qq annihilation 
[diagram (c) of Fig. [1] and quark-gluon scatter- 
ing, and the luminosity of quark-gluon states falls 
off more rapidly with the parton momenta than 
the luminosity of qq states [1, Q. The diphoton 



ME contributes 56% to the cross section, the pro- 
cesses qq — >■ 577 isr and gq 977isr 29%, and 
the process gq — >■ 977pgf^ 15%. Double radiation 
processes in minimum bias dijet events, such as 

99 ^ 997isr/fsr ^ISR/FSR ’ qq y 9*?Tisr/fsr ^ISR/FSR ’ 

99 ^ 997isr/fsr7isr/fsr ’ 99 ^ 997isr7isr ^nd 

gg — >■ 997fsr7fsrj were also examined but their 
overall contribution was estimated to only ^ 3% of 
the total, having no significant effect to any kine- 
matical distribution. Therefore, these processes 
were not included in the pythia calculation. 

• A fixed Next-to-Leading-Order (NLO) calculation 

using the diphox program Q. This generator ex- 
plictly includes parton fragmentation into photons 
[9|, i.e. processes in which nearly all the energy of 
a parton is transformed into a photon. Direct pro- 
duction contributes 85% to the cross section and 
fragmentation 15%. The diphox matrix element 
accounts for the qq — )■ 77 and gq — )■ 775 processes 
up to NLO, and LO for the gg 77 process, since 
this is already a second order process in the strong 
coupling. The NLO gg — ^ 77 contributions were 
examined with the GAMMA2mc program 0. Fig. 
HD shows an example of the uncorrected and cor- 
rected diphox predictions in comparison with the 
measured cross section as a function of the dipho- 
ton invariant mass M. The corrected prediction 
is calculated by running diphox without the LO 
gg 77 term and then adding the full LO-I-NLO 
99 77 calculation from GAMMA2mc incoher- 

ently, since the initial state is different in gluon fu- 
sion than in the other processes. The correction of 
the total cross section for the NLO gg — )■ 77 contri- 
bution is nearly 10%, which is comparable with the 
experimental and theoretical uncertainties (see Ta- 
ble [IV|. Therefore, this correction was not applied 
to the DIPHOX calculation. 

• A resummed NLO calculation using the RESBOS 
program 0. Here the effects of soft gluon ISR in 
the NLO calculation are analytically resummed to 
all orders in the strong coupling and reach NNLL 
accuracy. The resultant prediction is smoothly 
matched to the fixed-order NLO result in the kine- 
matic regions where the NLO matrix element is 
dominant. The RESBOS matrix element includes 
the qq — >■ 77, gq — )■ 779 and gg — >■ 77 processes up 
to NLO and it is adjusted so as to approximately 
account for fragmentation. 

All calculations are done by Monte Carlo event genera- 
tion and are subject to the experimental kinematic and 
isolation cuts. In the fixed-order NLO calculations the 
isolation cut is applied on parton variables and thus it 
only approximates the isolation cut appplied in the data 
and in pythia. The RESBOS predictions are restricted 
in the diphoton invariant mass M range from 2 rrib = 9 
GeV / to 2 mt = 350 GeV and they are shown up to 



16 



M = 300 GeV/c^ in the plots of the mass distribution, 
where mb and are the masses of the bottom and top 
quarks, respectively. 




default renormalization, factorization and fragmentation 
scales are all set to /r = M/2. In RESBOS the default 
renormalization and factorization scales are both set to 
/j, = M. In either case, all scales are varied by a fac- 
tor of 2 up and down relative to the default choice and 
this is taken as a conservative estimate of the total scale 
uncertainty. The proton PDF set is the CTEQ6.1M set 
for both DIPHOX and RESBOS. The corresponding 
uncertainty is estimated by varying the generated event 
weights within the 90% level uncertainties given by the 
20 CTEQ6.IM eigenvectors. 

The measured total cross section is shown in Table IIVI 
together with the predictions from the three theoretical 
calculations. All calculations are consistent with the size 
of the measured cross section within the experimental 
uncertainties. 



B. Kinematic Variables 




The complete description of the reaction hi+h 2 — >■ 71 -b 
72 -b A, where hi ^2 are hadrons, requires five independent 
kinematic variables. A suitable choice consists of the 
invariant mass 



M = ■y/2pT7lPT72[cOsh(y^i - 2/^2) - C0S((/)^1 - (j)j2)] 

( 2 ) 

the transverse momentum 

Pt = \Jpti 1 + PT72 + 2PT71PT72 cos((/)t,i - 4 >j 2 ) ( 3 ) 




FIG. 11 : The individual contributions to the cross section 
from events where both photons are generated according to 
the PYTHIA diphoton matrix element and from events where 
one photon originates from initial or final state radiation, as 
functions of the diphoton mass (top), transverse momentum 
(middle) and azimuthal difference (bottom). 



the rapidity 



PT7I coshy.yi -b PT72 coshy.y2 
the azimuthal difference 



A(/ = — (p^ 2 \ mod 7T (5) 

of the photon pair in the laboratory frame [3 , and the 
cosine of the pola r angle 9 of the 1®* photon in the Collins- 
Soper frame . This is defined as the rest frame of the 
photon pair chosen so that (a) the 3- momenta p/j and ph 2 
of the initial hadrons lie in the Oxz plane (with positive 
x) and (b) the z axis bisects the angle between^j,^ and 
—Ph 2 ■ This variable is generally determined by [33| 



2PT71PT72 sinh(y.^i - y^ 2 ) 

cos 0 = = — 

My/M^ + 



(6) 



NLO theoretical uncertainties are estimated for the 
choice of scale, representing the sensitivity to missing 
higher order terms, and for the PDEs. In diphox the 



For photons emitted at large angles with respect to the 
beam, cos 9 « tanh[(y.yi — 1 /^ 2 ) /2] in the limit Pt — >■ 0. In 
the above equations PT7i, y~/i and (f>yi are the transverse 
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FIG. 12: The measured cross section as a function of the 
diphoton mass in comparison with the diphox predictions 
without and with the NLO gg — >■ 77 correction, calculated 
by the GAMMA2 mc program. The shaded area is the total 
systematic uncertainty in the data. 

momentum, rapidity and azimuth of photon i, respec- 
tively, with i=l, 2. 

The set of {M, describes the kinematics of 

the diphoton system and, therefore, of possible heavy 
particles decaying into a photon pair, such as a Higgs 
boson. The existence of such a particle would mani- 
fest as a peak in the distribution of the invariant mass 
M. The results of this analysis are presented in the 
form of cross sections differential in each of the five kine- 
matic variables {M, Pr, A(f>,Yyj,cos9} and in the vari- 
able z = PT72/PT7I) the ratio of sub-leading to lead- 
ing photon transverse momentum (0<z<l). Three kine- 
matic cases are examined: 

• Differential cross sections without additional kine- 
matic cuts. No kinematic cut other than those 
listed in Table imi is applied. The results of this 
case are presented in Subsection IV Cl 

• Differential cross sections for Pt<M. The kine- 
matics in this case are similar to the diphoton de- 
cay of a heavy particle, such as a Higgs boson, pro- 
duced in events of moderate parton activity. At the 
Tevatron, prompt photon pairs are almost entirely 
produced in this case by 1ow-Pt quark-antiquark 
annihilation. The results of this case are presented 
in Subsection IV PI 

• Differential cross sections for Pt>M. The impor- 
tance of high-Pr contributions from gluon-gluon fu- 
sion, fragmentations and ISR is enhanced in this 
case. The results of this case are presented in Sub- 
section IV El 

The oveflow data entries are excluded from the M and 
Pt histograms, to keep the cross section definition con- 
sistent for the data and the theories at the highest bins. 



TABLE IV: The total diphoton production cross section ob- 
tained from the measurement and from the theoretical calcu- 
lations. The PYTHIA 77 calculation involves only the qq — > 77 
and gg — >■ 77 processes. The pythia 77 -|- 7] calculation in- 
cludes also the qq — >■ 773 and gq — >■ 775 processes. 



Cross section (pb) 


Data 


12.47 ± 0.2Etat ±3.74yst 


RESBOS 


11.31 ± 2.45syst 


DIPHOX 


10.58 ± 0.55syst 


PYTHIA 77 -1- 7j 


9.19 


PYTHIA 77 


5.03 



For each kinematic variable, the following plots are pre- 
sented: 

• The measured and calculated cross sections as func- 
tions of the selected variable. Each of these plots 
includes the predictions of all three calculations, 
for comparison, and shows only the uncertainties 
of the data. The prediction of the pythia 77 cal- 
culation, involving only the qq —I 77 and gg —I 77 
processes, is also shown in these plots, to be com- 
pared with the PYTHIA 77-l-7j calculation (j=jet), 
which includes also the qq — i 773 and gq — i 779 
processes. 

• The relative deviations of the data from each calcu- 
lation, in the form (data— theory) /theory, as func- 
tions of the selected variable. These plots show the 
comparison of the data with each calculation sep- 
arately and include the uncertainties of the NLO 
predictions. No relative deviations are shown for 
the PYTHIA 77 calculation. The benchmark parton 
showering MC calculation, compared in detail with 
the data, is pythia 77-1-7). 

C. Differential Cross Sections Without Additional 
Kinematic Cut 

Fig. IT^ shows the results for da/dM and da/dP^. The 
mass spectrum peaks at M = 2 a / ”«32 GeV/c^. 
All three theoretical predictions for da/dM are in rea- 
sonable agreement with the data, within uncertainties, 
except in the region 6 GeV/c^<M <32 GeV/c^ below the 
mass peak. The low mass limit of 6 GeV /(? is set by the 
Ai?(77) and cuts. This region is rich in events com- 
ing from gluon scattering and fragmentation. All three 
predictions underestimate the data in this region. 

The excess of the data over all three predictions for M 
below the peak of the mass spectrum is reflected in the re- 
gion 20 GeV /c<Pt<50 GeV /c of the Pt spectrum, which 
has a shoulder around Pt=Pt 2 ”+Pti "=32 GeV/c (the 
so-called “Guillet shoulder” ) . This arises from a collinear 
enhancement for the two photons in the fragmentation 
processes which, however, is suppressed by the 
cut. The RESBOS predictions for da/dP^ are in overall 
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FIG. 13: The cross section as a function of the diphoton invariant mass (left) and transverse momentum (right). Top: the 
absolute cross section values. Bottom: the relative deviations of the data from the predictions. Note: the vertical axes scales 
differ between relative deviation plots. The shaded area is the total systematic uncertainty in the data. 



agreement with the data, within uncertainties, except in 
this region. The diphox prediction underestimates the 
data, in addition, for Pt< 20 GeV/c, where the resum- 
mation effects implemented in resbos provide a better 
description. The pythia prediction underestimates the 
data at very low Pt, Pt< 10 GeV /c, showing that the LL 
resummation of parton showering is less accurate than 
the NNLL resummation implemented in RESBOS. The 
PYTHIA prediction is in reasonably good agreement with 
the data, within uncertainties, in the rest of the Pp range 
due to the qq —I 577 isji and gq —I q^%g-n processes which 
make the pythia Pp spectrum sufficiently hard (see the 
middle plot of Fig. ITT]) . 

Fig. [HI shows the results for da/dA(f> and da/dY^^. 



The A(j) spectrum peaks at A(^= 7 t, corresponding to van- 
ishing diphoton Pp, and the spectrum at Fy-,,=0, 
corresponding to vanishing diphoton momentum P^ par- 
allel to the proton beam. While all three predictions 
agree fairly well with the measured da/dY^-y, within un- 
certainties, all three of them underestimate the data in 
the low end of the A(j) spectrum. This region is domi- 
nated by events with low mass and high Pp . The RESBOS 
prediction provides the best description of the measured 
da / dA(j) for A(p>2.5 rad, where soft gluon resummation 
is important, pythia provides the best description in 
the region 1 rad<A(/)<2.5 rad where the qq —I 577isr 
and gq —I 97713 ^ processes are the most important (see 
the bottom plot of Fig. HH). 
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FIG. 14: The cross section as a function of the diphoton azimuthal distance (left) and of the diphoton rapidity (right). Top: 
the absolute cross section values. Bottom: the relative deviations of the data from the predictions. Note: the vertical axes 
scales differ between relative deviation plots. The shaded area is the total systematic uncertainty in the data. 



Fig. [T5] shows the results for da/dcosO and da/dz. All 
three predictions agree with the data, within uncertain- 
ties. Exceptions are the predictions of all three calcu- 
lations underestimating the data in the two ends of the 
cos9 spectrum, where again gluon scattering processes 
and associated fragmentation are expected to dominate 

i- 

In general, all three calculations reproduce most of 
the main features of the data, as observed in the ear- 
lier diphoton cross section measurements Eini- How- 
ever, depending on their approximations, they display 
differences with each other and with the data in certain 
kinematic regions. There is a problem common to all 
three calculations in the description of events with very 



low diphoton mass, low azimuthal distance and diphoton 
transverse momentum in the region of the “Guillet shoul- 
der” . Such events include fragmentation at a relatively 
high rate. The pythia 77 calculation fails completely to 
describe the data both in the scale, where it is low by 
a factor of 2.5, and in the shape, particularly of the Pt, 
A(^ and z distributions, where it predicts a much softer 
spectrum than the data. This is in agreement with the 
conclusion of Ref. [l^ which tested only pythia 77 as 
a parton showering MC prediction. 
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FIG. 15: The cross section as a function of the cosine of the polar angle in the Collins-Soper frame (left) and of the ratio of the 
subleading photon Et to leading photon Et (right). Top: the absolute cross section values. Bottom: the relative deviations 
of the data from the predictions. Note: the vertical axes scales differ between relative deviation plots. The shaded area is the 
total systematic uncertainty in the data. 



D. Differential Cross Sections for Pt<M 
Kinematics 

Fig. ITBl shows the results for da/dM and da/dPr for 
Pt<M. The low tail of the mass spectrum, in the region 
6 GeV/c^<M<32 GeV/c^, and the shoulder of the Pt 
spectrum, in the region 20 GeV/c<PT<50 GeV/c, are 
now eliminated. The agreement between the data and 
all three predictions is improved in this case. However, 
DIPHOX still underestimates the data for Pt<40 GeV/c 
and similarly pythia still underestimates the data for 
Ft <10 GeV/c, thus showing the importance of NNLL 
low- F t resummation in this case as well. 



Fig. [T 7 I shows the results for da/dAcj) and da/dY^^ for 
Pt:<M . The tail of the A</) spectrum for A^<7 t/ 2 is now 
weaker but the measured cross section is underestimated 
by all three predictions, as in the case of unconstrained 
kinematics. 

Fig. [TSl shows the results for da /d cos 9 and da/dz 
for Pt<M. The results are similar to the case of un- 
constrained kinematics. Generally, all three calculations 
agree with the data, within uncertainties. Exceptions are 
again the predictions of all three calculations in the two 
ends of the cos 9 spectrum, where they underestimate the 
data. 

In general, events with kinematics similar to the de- 
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FIG. 16: The cross section as a function of the diphoton invariant mass (left) and transverse momentum (right) for Pt<M. 
Top: the absolute cross section values. Bottom: the relative deviations of the data from the predictions. Note: the vertical 
axes scales differ between relative deviation plots. The shaded area is the total systematic uncertainty in the data. 



cay of a heavy particle with low transverse momentum 
into a photon pair, such as gg— i-H— >77 production and 
decay, are better described by the theory than events 
with low mass and hi gh t ransverse momentum. This is 
also observed in Ref. [ij which examines only the case 
of Pt<M. This observation is important for current 
searches of yet undiscovered particles with a diphoton 
decay signature. The pythia 77 calculation again fails 
to describe the data both in the scale and in the shape, in 
agreement with the conclusion of Ref. m which tested 
only PYTHIA 77 as a parton showering MC prediction. 



E. Differential Cross Sections for Pt>M 
Kinematics 

Fig. [T^ shows the results for da/dM and da/dP^ for 
Pt'>M. Both spectra are depleted in this case: the mass 
spectrum for M>200 GeV/c^ and the transverse momen- 
tum spectrum for Pt< 20 GeV/c. All three calculations 
underestimate the data. 

Fig. [20I shows the results for da/dAcj) and dajdY^^ 
for Pt>M. The A(f> spectrum is strongly suppressed for 
A(^>7 t/ 2. Again, the measured cross section is underes- 
timated by all three calculations. 

Fig. shows the results for da /d cos 9 and da/dz for 
Pt>M. In contrast with the unconstrained kinematics 
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FIG. 17: The cross section as a function of the diphoton azimuthal distance (left) and of the diphoton rapidity (right) for 
Pt<M. Top: the absolute cross section values. Bottom: the relative deviations of the data from the predictions. Note: the 
vertical axes scales differ between relative deviation plots. The shaded area is the total systematic uncertainty in the data. 



and the Pt<M kinematics, in this case all three calcula- 
tions underestimate the data through the full ranges of 
the cos 9 and z spectra. 



In general, events with low diphoton mass and high 
diphoton transverse momentum, mainly coming from 
fragmentation, are not well described by the examined 
calculations. This observation is important for measure- 
ments under conditions where contributions from such 
events are strong, as in the LHC [l|. 



VI. CONCLUSIONS 

In summary, the prompt diphoton production cross 
section, differential in kinematic variables sensitive to the 
dynamics of the reaction mechanism, is measured using 
data corresponding to an integrated of luminosity 5.36 
fb“^ collected with the CDF II detector. The large size 
of the data sample allows for scanning a much more ex- 
tended phase space and with a better statistical precision 
than in earlier measurements. Using a novel technique 
for the background subtraction, based on the track iso- 
lation, the overall systematic uncertainty is limited to 
about 30% on average. 

The results of the measurement are compared with 
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FIG. 18: The cross section as a function of the cosine of the polar angle in the Collins-Soper frame (left) and of the ratio of the 
subleading photon Et to leading photon Et (right) for Pt<M. Top: the absolute cross section values. Bottom: the relative 
deviations of the data from the predictions. Note: the vertical axes scales differ between relative deviation plots. The shaded 
area is the total systematic uncertainty in the data. 



three state-of-the-art calculations, applying complemen- 
tary techniques in modeling the reaction. All three cal- 
culations describe events with large diphoton mass and 
small diphoton transverse momentum fairly well, where 
the kinematics is similar to the decay of a 1ow-Pt heavy 
particle, such as the Higgs boson, decaying into a photon 
pair. Exceptions are kinematic regions where gluon inter- 
actions and the associated fragmentations of quarks into 
photons are expected to be important, such as the low 
mass and azimuthal difference regions and the region of 
the “Guillet shoulder” at moderate transverse momen- 
tum. All three calculations underestimate the data in 
those regions. Although the diphox calculation explic- 



itly includes a fragmentation model, it fails to reproduce 
the data in those sensitive regions, possibly because of 
the approximate nature of the requirement of photon iso- 
lation in the diphox framework. This requirement is 
mostly responsible for the suppression of fragmentation 
contributions and is applied using hadron variables in 
the data but using parton variables in diphox. The low 
transverse momentum and large azimuthal difference re- 
gions, where resummation in the diphoton transverse mo- 
mentum is important, are best described by RESBOS, as 
expected from the analytical resummation implemented 
in this calculation. Photon radiation, especially from the 
initial state quarks, in addition to the prompt photon 
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FIG. 19: The cross section as a function of the diphoton invariant mass (left) and transverse momentum (right) for Pt>M. 
Top: the absolute cross section values. Bottom: the relative deviations of the data from the predictions. Note: the vertical 
axes scales differ between relative deviation plots. The shaded area is the total systematic uncertainty in the data. 



production at the hard scattering, is for the first time 
shown to play a very important role in the parton show- 
ering PYTHIA calculation in order to bring the prediction 
into reasonable agreement with the data. 
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FIG. 20: The cross section as a function of the diphoton azimuthal distance (left) and of the diphoton rapidity (right) for 
Pt>M. Top: the absolute cross section values. Bottom: the relative deviations of the data from the predictions. Note: the 
vertical axes scales differ between relative deviation plots. The shaded area is the total systematic uncertainty in the data. 
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Appendix A: Content likelihood of the data sample 

The signal and background disentanglement is done in 
a maximum likelihood framework. A likelihood of the 
composition of the baseline 77 sample is defined on the 
basis of binomial probabilities for the observed photons 
to pass or fail the track isolation cut: 

^ = (AI) 

ik 



where fc = I, . . . , A labels all events in the sample and i = 
PP-: pf-: fPy ff labels the categories of events in which both 
photons pass, the leading passes and the sub-leading fails, 
the leading fails and the sub-leading passes, and both 
photons fail the track isolation cut, respectively, /ii(fc) 
is the probability of the event k to fall in the category i 
and rii{k) is the observation for the event k; i.e. rii(k) 
is 1 for one of the four categories i and 0 for the other 
three. The probabilities can be analyzed as follows: 

P^{k)=Y^ Eij {k)pj (fc) ( A2) 

3 

where j = ss, sb, bs, bb labels the categories in which both 
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FIG. 21: The cross section as a function of the cosine of the polar angle in the Collins-Soper frame (left) and of the ratio of the 
subleading photon Et to leading photon Et (right) for Pt>M. Top: the absolute cross section values. Bottom: the relative 
deviations of the data from the predictions. Note: the vertical axes scales differ between relative deviation plots. The shaded 
area is the total systematic uncertainty in the data. 



photons are signal, the leading is signal and the sub- 
leading background, the leading is background and the 
sub-leading signal, and both photons are background, re- 
spectively. Pj{k) is the probability for an event k to be 
in the category j. Sij{k) is the probability for an event k 
of the truth category j to be observed in the observation 



category i. This probability is directly related with the 
efficiencies of the two photons in the event k to pass the 
track isolation cut, since the efficiencies are defined as 
the probabilities for the leading or sub-leading photon to 
pass the cut. In a 4 x 4 matrix notation. 



e{k) 



esi{k)€s2{k) 

€sl{k){l - €s2{k)) 

(I - esi{k))€s2{k) 

(1 - esi(fc))(I - es 2 (fc)) (1 



esi{k)eb2{k) 

€siik){l - eb2{k)) 

(I - €sl{k))€b2{k) 

- esi(fc))(I - eb 2 (fc)) (1 



ebi{k)es2{k) 

£bl(fc)(l - €s2{k)) 

(1 - €bl{k))€s2{k) 

- ebi(fc))(I - es 2 (fc)) (1 



ebi{k)eb2{k) 
£bl(fc)(l - £b2(fc)) 

(1 - £bl(fc))£&2(fc) 

- £bl(fc))(I - £62(fc)) 



(A3) 
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where eap{k) (a = s or 6, /3 = 1 or 2) are the efficiencies 
of the leading or sub-leading photon, coming from the 
signal or from the background, to pass the track isolation 
cut. 

The likelihood C is maximized or, equivalently, the op- 
posite of its natural logarithm 



L = -ln£ (A4) 

= ^ {ni{k) In^i(fc) -I- [1 - ni{k)] In [1 - ^J,i{k)] } 

ik 

is minimized with respect to the probabilities Pj(k). The 
minimization of the logarithm L leads to the system of 
equations 



m{k) = pi{k) ='^ Eij {k)pj (k) ( A5) 

3 

with solutions 

P] W ik)n^ (k) ( A6) 

i 

The choice of the track isolation cut at 1 GeV / c gives the 
efficiencies Capik) sufficient discriminating power among 
the truth categories j for the matrix e{k) to be non- 
singular. By summing the probabilties over all events 
in the baseline sample the maximum likelihood composi- 
tion of the sample is obtained: 



Wj = '^Pj{k) = (A7) 

k ik 

Eq. IA7I for the truth category j = ss provides the signal 
fraction in the baseline 77 sample. 

In general, the composition of a sample of events with 
m photons each can be resolved in a maximum likeli- 
hood framework by the inversion of a 2™ x 2™ matrix, 
constructed as in Eq. IA3i which transforms the prob- 
ability 2’”-vectors Pj{k) to the observation 2’”-vectors 
rii{k). The generic matrix element £ij{k) contains a fac- 
tor €ap{k) or 1 — Capik) for each photon /3 = 1, 2, . . . , m 
with a = s {b) if the photon is signal (background) in 



the truth category j and passes or fails, respectively, the 
track isolation cut in the observation category i. In this 
context, Eq. [T] for the single photon sample is derived 
from Eq. IA7I by inverting for each event the matrix 



£T = 



is Cb \ 

1 - Cs 1- €b J 



=> £ 



-1 



1 ( 1-et 

es- £b V 




with observation vectors 



rip = 





(A8) 



(A9) 



for photons passing (rip) or failing (n/) the track isolation 
cut. Eq. IA7l then gives for the signal fraction of the single 
photon sample 




where, in the first line, the first sum runs over all photons 
passing the cut and the second sum runs over all photons 
failing the cut. Eq. lAlOl is identical with Eq. [T] 



Appendix B: Cross Section Tables 

This Appendix provides a set of Tables of the measured 
differential cross section values as functions of the six 
kinematic variables selected in this analysis. Each Table 
lists the bins of the selected variable, the values of the 
cross section in the respective bins for the three examined 
cases of no kinematic cut, Pt<M, and Pt>M, and the 
statistical and total systematic uncertainties associated 
with each cross section value. 
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TABLE V: The diphoton production cross section differential in the diphoton invariant mass. The first error in the cross section 



is statistical and the second systematic. 


Mass bin [GeV/c'^] 


Cross section without cut 


Cross section for Pt<M 


Cross section for Pt>M 




]pb/(GeV/c^)] 


[pb/(GeV/c^)] 




[pb/(GeV/c^)] 


0 - 5 






0 






0 






0 






5 - 10 


0.014679 


± 


0.002997 


± 


0.003369 


0 


0.014679 


± 


0.002997 


± 


0.003301 


10 - 15 


0.004441 


± 


0.000607 


± 


0.001115 


0 


0.036310 


± 


0.004965 


± 


0.009118 


15 - 20 


0.027568 


± 


0.005232 


± 


0.010986 


0 


0.029520 


± 


0.005602 


± 


0.011764 


20 - 25 
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± 
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± 
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± 
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± 
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± 
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± 
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± 
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± 
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± 
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± 
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± 
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± 
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± 
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± 
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40 - 45 
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± 
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45 - 50 
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0.088299 
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0.003078 


± 
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± 
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50 - 55 
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0.011824 


± 
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0.201493 ± 0.011730 ± 0.056897 


0.005372 


± 


0.001299 


± 
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± 
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± 


0.036287 
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0.003014 


± 
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± 
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60 - 65 
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± 
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0.001071 


± 
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± 
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± 
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± 
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± 
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± 
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± 
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± 
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± 
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± 
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± 
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0.000480 
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± 
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± 
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± 
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TABLE VI: The diphoton production cross section differential in the diphoton transverse momentum. The first error in the 
cross section is statistical and the second systematic. 



Pt bin ]GeV/c] 


Cross section without cut 


Cross section for Pt<M 


Cross section for Pt>M 




[pb/(GeV/c)] 


]pb/(GeV/c)] 


[pb/(GeV/c)] 



0 - 1 


0.276549 


± 


0.030968 


± 


0.067613 


0.276542 


± 


0.030967 


± 


0.067612 


0 






1 - 2 


0.823081 


± 


0.051902 


± 


0.184444 


0.823132 


± 


0.051905 


± 


0.184456 


0 






2 - 3 


1.043320 


± 


0.056713 


± 


0.236673 


1.043373 


± 


0.056716 


± 


0.236685 


0 






3 - 4 


1.026369 


± 


0.055895 


± 


0.241145 


1.026574 


± 


0.055907 


± 


0.241193 


0 






4 - 5 


0.913944 


± 


0.055127 


± 


0.234272 


0.914185 


± 


0.055142 


± 


0.234334 


0 






5 - 6 


0.908628 


± 


0.053530 


± 


0.215222 


0.909148 


± 


0.053560 


± 


0.215345 


0 






6 - 7 


0.784121 


± 


0.049343 


± 


0.200779 


0.784624 


± 


0.049375 


± 


0.200908 


0 






7 - 8 


0.606111 


± 


0.045730 


± 


0.166151 


0.606435 


± 


0.045755 


± 


0.166240 


0 






8 - 9 


0.538172 


± 


0.042695 


± 


0.159898 


0.538651 


± 


0.042733 


± 


0.160040 


0 






9 - 10 


0.416951 


± 


0.039170 


± 


0.125580 


0.417071 


± 


0.039181 


± 


0.125616 


0 






10 - 12 


0.371872 


± 


0.025701 


± 


0.113307 


0.372199 


± 


0.025723 


± 


0.113406 


0 






12 - 14 


0.275282 


± 


0.022571 


± 


0.089464 


0.275797 


± 


0.022613 


± 


0.089631 


0 






14 - 16 


0.196349 


± 


0.019572 


± 


0.072611 


0.196591 


± 


0.019597 


± 


0.072700 


0 






16 - 18 


0.179301 


± 


0.017290 


± 


0.060538 


0.179549 


± 


0.017314 


± 


0.060622 


0 






18 - 20 


0.125584 


± 


0.015368 


± 


0.045351 


0.126581 


± 


0.015490 


± 


0.045711 


0 






20 - 25 


0.127625 


± 


0.009118 


± 


0.038128 


0.125027 


± 


0.009019 


± 


0.037597 


0.001932 ± 0.000883 


± 


0.000457 


25 - 30 


0.092613 


± 


0.007839 


± 


0.028973 


0.078351 


± 


0.006931 


± 


0.024441 


0.013751 ± 0.003737 


± 


0.004417 


30 - 35 


0.087187 


± 


0.007329 


± 


0.027382 


0.040507 


± 


0.004874 


± 


0.012880 


0.049065 ± 0.005765 


± 


0.015309 


35 - 40 


0.056515 


± 


0.006341 


± 


0.020793 


0.024564 


± 


0.003619 


± 


0.007531 


0.033546 ± 0.005550 


± 


0.014138 


40 - 45 


0.039824 


± 


0.004995 


± 


0.012246 


0.012094 


± 


0.002774 


± 


0.004861 


0.030706 ± 0.004582 


± 


0.008315 


45 - 50 


0.022650 


± 


0.003604 


± 


0.006110 


0.010569 


± 


0.002027 


± 


0.002311 


0.012983 ± 0.003290 


± 


0.004252 


50 - 60 


0.018204 


± 


0.002042 


± 


0.004416 


0.006578 


± 


0.001173 


± 


0.001669 


0.012472 ± 0.001798 


± 


0.002971 


60 - 70 


0.007542 


± 


0.001410 


± 


0.002019 


0.003080 


± 


0.000766 


± 


0.000693 


0.004862 ± 0.001318 


± 


0.001481 


70 - 80 


0.005717 


± 


0.001043 


± 


0.001271 


0.001773 


± 


0.000483 


± 


0.000331 


0.004388 ± 0.001044 


± 


0.001080 


80 - 90 


0.003467 


± 


0.000656 


± 


0.000676 


0.000801 


± 


0.000292 


± 


0.000167 


0.003387 ± 0.000749 


± 


0.000642 


90 - 100 


0.001355 


± 


0.000487 


± 


0.000377 


0.000079 


± 


0.000188 


± 


0.000103 


0.001242 ± 0.000433 


± 


0.000279 


100 - 120 


0.001057 


± 


0.000212 


± 


0.000281 


0.000198 


± 


0.000097 


± 


0.000062 


0.000972 ± 0.000212 


± 


0.000239 


120 - 140 


0.000655 


± 


0.000228 


± 


0.000135 


0.000059 


± 


0.000071 


± 


0.000017 


0.000516 ± 0.000188 


± 


0.000094 


140 - 160 


0.000287 


± 


0.000122 


± 


0.000051 


0.000008 


± 


0.000030 


± 


0.000007 


0.001222 ± 0.000496 


± 


0.000227 


160 - 200 


0.000100 


± 


0.000047 


± 


0.000021 






0 






0.000227 ± 0.000105 


± 


0.000039 



~87% of the total area. 

[24] T. Aaltonen et al. (CDF Collaboration), Phys. Rev. D 
82, 052005 (2010). 

[25] E. Gerchtein and M. Paulini, CHEP-2003-TUMT005, 
arXiv:physics/0306031 (2003); the version of GEANT used 
for the detector simulation is 3.21, see the CERN Pro- 
gram Library Long Writeup W5013. 

[26] T. Aaltonen et al. (CDF Collaboration), Phys. Rev D 77, 
092001 (2008); Phys. Rev. Lett. 102, 232002 (2009). 

[27] H. L. Lai, J. Huston, S. Kuhlmann, J. Morhn, F. Olness, 
J.F. Owens, J. Pumplin and W. K. Tung, Eur. Phys. J. 
C 12, 375 (2000). 



[28] R. Field and C. Group, arXiv:0510198 ]hep-ph] (2005). 

[29] T. Aaltonen et al. (CDF Collaboration), Phys. Rev. D 
80, 111106 (2009). 

[30] W. M. Yao et al., “Review of Particle Physics”, J. Phys. 
G 33, 2006. 

[31] D. Acosta et al. (GDF Collaboration), Phys. Rev. Lett. 
94, 091803 (2005). 

[32] J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P. Nadol- 
sky and W. K. Tung, J. HEP 0207, 012 (2002). 

[33] J. C. Collins and D. E. Soper, Phys. Rev D 16, 2219 
(1977). 



30 



TABLE VII: The diphoton production cross section differential in the diphoton azimuthal difference. The first error in the 
cross section is statistical and the second systematic. 

A(j) bin [radians] Cross section without cut Cross section for Pt<M Cross section for Pt>M 

[pb/rad] [pb/rad] [pb/rad] 

0.000 - 0.105 0.456343 ± 0.099393 ± 0.099826 0.008959 ± 0.006342 ± 0.004135 0.455591 ± 0.101696 ± 0.101433 

0.105 - 0.209 0.443854 ± 0.103665 ± 0.101738 0 0.506682 ± 0.116125 ± 0.114737 

0.209 - 0.314 0.346811 ± 0.090621 ± 0.085583 0.000964 ± 0.000683 ± 0.000446 0.377053 ± 0.101504 ± 0.096562 

0.314 - 0.419 0.805494 ± 0.112525 ± 0.174161 0.008695 ± 0.006135 ± 0.001692 0.835398 ± 0.118774 ± 0.181144 

0.419 - 0.524 0.621352 ± 0.134650 ± 0.162677 0 0.668658 ± 0.139966 ± 0.168432 

0.524 - 0.628 0.695012 ± 0.138448 ± 0.183259 0.007809 ± 0.011151 ± 0.001958 0.712287 ± 0.143289 ± 0.188042 

0.628 - 0.733 0.553260 ± 0.152492 ± 0.264377 0.023257 ± 0.020326 ± 0.006509 0.528230 ± 0.151279 ± 0.258421 

0.733 - 0.838 0.375583 ± 0.144617 ± 0.261528 0 0.433524 ± 0.158014 ± 0.270138 

0.838 - 0.942 0.671584 ± 0.143112 ± 0.228641 0.029115 ± 0.020067 ± 0.010039 0.679242 ± 0.151863 ± 0.231038 

0.942 - 1.047 0.522112 ± 0.137889 ± 0.219688 0.038178 ± 0.018111 ± 0.008311 0.496178 ± 0.145545 ± 0.230738 

1.047 - 1.152 0.864793 ± 0.144971 ± 0.220201 0.087654 ± 0.032029 ± 0.018939 0.781381 ± 0.145224 ± 0.207021 

1.152 - 1.257 0.798674 ± 0.154962 ± 0.224791 0.150640 ± 0.043789 ± 0.033689 0.638090 ± 0.150851 ± 0.210793 

1.257 - 1.361 1.124851 ± 0.168956 ± 0.257447 0.188806 ± 0.056573 ± 0.044931 0.939944 ± 0.162249 ± 0.212925 

1.361 - 1.466 0.639906 ± 0.169993 ± 0.284205 0.118937 ± 0.079298 ± 0.074899 0.527203 ± 0.151112 ± 0.210203 

1.466 - 1.571 0.834419 ± 0.171606 ± 0.270175 0.421297 ± 0.113865 ± 0.107330 0.408768 ± 0.128128 ± 0.171398 

1.571 - 1.676 1.250665 ± 0.169418 ± 0.311391 1.007357 ± 0.159164 ± 0.278003 0.230173 ± 0.057530 ± 0.047005 

1.676 - 1.780 0.832275 ± 0.186990 ± 0.364452 0.826622 ± 0.181796 ± 0.334001 0.006883 ± 0.044420 ± 0.000755 

1.780 - 1.885 1.001287 ± 0.187680 ± 0.340539 0.998464 ± 0.185177 ± 0.326921 0.004982 ± 0.030035 ± 0.000548 

1.885 - 1.990 1.076470 ± 0.209360 ± 0.442750 0.985068 ± 0.206856 ± 0.435764 0.097539 ± 0.033960 ± 0.017947 

1.990 - 2.094 1.667315 ± 0.233681 ± 0.501973 1.628501 ± 0.232554 ± 0.498599 0.030895 ± 0.019985 ± 0.005693 

2.094 - 2.199 2.275603 ± 0.257154 ± 0.585295 2.271233 ± 0.256406 ± 0.580308 0.005091 ± 0.018432 ± 0.007031 

2.199 - 2.304 1.982347 ± 0.275434 ± 0.703905 1.945552 ± 0.274931 ± 0.703233 0.021202 ± 0.011392 ± 0.004164 

2.304 - 2.409 3.180244 ± 0.298472 ± 0.883152 3.185517 ± 0.298001 ± 0.878647 0 

2.409 - 2.513 3.344205 ± 0.354571 ± 1.099021 3.345126 ± 0.354498 ± 1.097423 0.000142 ± 0.006657 ± 0.000016 

2.513 - 2.618 4.913562 ± 0.403046 ± 1.447020 4.883607 ± 0.402930 ± 1.445024 0.010602 ± 0.005301 ± 0.001996 

2.618 - 2.723 6.787434 ± 0.476941 ± 1.881799 6.788491 ± 0.477385 ± 1.881996 0.000988 ± 0.001690 ± 0.000302 

2.723 - 2.827 9.949192 ± 0.584235 ± 2.773900 9.948025 ± 0.584526 ± 2.773673 0.001445 ± 0.002278 ± 0.000396 

2.827 - 2.932 14.781949 ± 0.694936 ± 3.966263 14.791632 ± 0.695391 ± 3.968861 0 

2.932 - 3.037 21.597660 ± 0.861897 ± 5.826928 21.602392 ± 0.862322 ± 5.828225 0.000218 ± 0.000352 ± 0.000061 

3.037 - 3.142 33.827076 ± 0.998489 ± 8.203793 33.825920 ± 0.998661 ± 8.204532 0.000133 ± 0.000133 ± 0.000024 
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TABLE VIII: The diphoton production cross section differential in the diphoton rapidity. The first error in the cross section is 
statistical and the second systematic. 

Ygg bin Cross section without cut Cross section for Pt<M Cross section for Pt>M 



[py by bM 



- 1.20 - - 1.12 


0.047796 


± 


0.032826 


± 


0.014355 


0.047253 


± 


0.027906 


± 


0.017471 


0.000543 


± 


0.017286 


± 


0.000059 


- 1.12 - - 1.04 


0.357094 


± 


0.120949 


± 


0.102530 


0.353126 


± 


0.113517 


± 


0.089113 


0.003968 


± 


0.041744 


± 


0.000434 


- 1.04 - - 0.96 


0.187965 


± 


0.037254 


± 


0.047292 


0.179976 


± 


0.034997 


± 


0.041721 


0.004941 


± 


0.006801 


± 


0.004111 


- 0.96 - - 0.88 


1.216128 


± 


0.179918 


± 


0.350105 


1.123226 


± 


0.169654 


± 


0.303864 


0.087621 


± 


0.056102 


± 


0.049798 


- 0.88 - - 0.80 


2.934750 


± 


0.300535 


± 


0.734394 


2.634538 


± 


0.284117 


± 


0.653983 


0.284037 


± 


0.092706 


± 


0.076637 


- 0.80 - - 0.72 


3.523381 


± 


0.364743 


± 


1.002344 


3.275366 


± 


0.344374 


± 


0.929432 


0.243582 


± 


0.121183 


± 


0.071961 


- 0.72 - - 0.64 


4.727782 


± 


0.415361 


± 


1.224628 


4.251420 


± 


0.388978 


± 


1.079371 


0.442221 


± 


0.136633 


± 


0.142010 


- 0.64 - - 0.56 


5.043031 


± 


0.463117 


± 


1.552137 


4.391841 


± 


0.437811 


± 


1.375512 


0.704898 


± 


0.160548 


± 


0.190566 


- 0.56 - - 0.48 


6.849294 


± 


0.521381 


± 


1.783643 


6.156165 


± 


0.496178 


± 


1.631273 


0.704821 


± 


0.161752 


± 


0.160419 


- 0.48 - - 0.40 


7.515287 


± 


0.565874 


± 


1.957807 


6.894050 


± 


0.542207 


± 


1.818195 


0.612712 


± 


0.159704 


± 


0.141450 


- 0.40 - - 0.32 


7.432176 


± 


0.605607 


± 


2.230531 


7.009432 


± 


0.582405 


± 


2.029948 


0.387476 


± 


0.149935 


± 


0.191711 


- 0.32 - - 0.24 


9.473290 


± 


0.644927 


± 


2.438148 


9.018377 


± 


0.624024 


± 


2.249574 


0.448143 


± 


0.156451 


± 


0.207981 


- 0.24 - - 0.16 


9.171682 


± 


0.650452 


± 


2.482947 


8.178396 


± 


0.625007 


± 


2.276978 


1.030157 


± 


0.186527 


± 


0.232884 


- 0.16 - - 0.08 


10.391827 ± 0.673024 ± 2.604154 


9.833859 


± 


0.649550 


± 


2.409783 


0.581114 


± 


0.183635 


± 


0.235125 


- 0.08 - 0 


10.685762 ± 0.683909 ± 2.779428 


9.737009 


± 


0.655963 


± 


2.536669 


1.011320 


± 


0.205173 


± 


0.258684 


0 - 0.08 


9.778326 ± 0.680032 ± 2.746355 


8.910962 


± 


0.656693 


± 


2.541139 


0.779653 


± 


0.160122 


± 


0.189297 


0.08 - 0.16 


10.245283 ± 0.664715 ± 2.572113 


9.428201 


± 


0.642795 


± 


2.382707 


0.844021 


± 


0.173754 


± 


0.196508 


0.16 - 0.24 


8.898950 


± 


0.650517 


± 


2.522506 


8.174833 


± 


0.625682 


± 


2.299986 


0.740208 


± 


0.181740 


± 


0.229029 


0.24 - 0.32 


9.873399 


± 


0.639878 


± 


2.378116 


8.881087 


± 


0.614339 


± 


2.143591 


0.980913 


± 


0.177119 


± 


0.232361 


0.32 - 0.40 


8.208516 


± 


0.599002 


± 


2.150160 


7.608601 


± 


0.571739 


± 


1.962087 


0.575259 


± 


0.170775 


± 


0.186692 


0.40 - 0.48 


8.310832 


± 


0.582531 


± 


2.114900 


7.805289 


± 


0.557393 


± 


1.950114 


0.489608 


± 


0.162967 


± 


0.168822 


0.48 - 0.56 


6.158322 


± 


0.525197 


± 


1.682342 


5.551052 


± 


0.501466 


± 


1.536281 


0.602316 


± 


0.155006 


± 


0.146833 


0.56 - 0.64 


5.764018 


± 


0.484449 


± 


1.615546 


5.254852 


± 


0.463164 


± 


1.480387 


0.515143 


± 


0.143704 


± 


0.136823 


0.64 - 0.72 


3.746496 


± 


0.413343 


± 


1.272077 


3.516958 


± 


0.394670 


± 


1.144314 


0.219857 


± 


0.116288 


± 


0.126414 


0.72 - 0.80 


2.638815 


± 


0.344297 


± 


0.937298 


2.365470 


± 


0.327398 


± 


0.833140 


0.236348 


± 


0.092090 


± 


0.090185 


0.80 - 0.88 


1.908901 


± 


0.278836 


± 


0.661709 


1.745827 


± 


0.263054 


± 


0.558774 


0.178584 


± 


0.101977 


± 


0.123069 


0.88 - 0.96 


1.342328 


± 


0.185391 


± 


0.348576 


1.092484 


± 


0.174767 


± 


0.311890 


0.266332 


± 


0.065150 


± 


0.055395 


0.96 - 1.04 


0.125836 


± 


0.034791 


± 


0.043927 


0.121207 


± 


0.033430 


± 


0.042727 


0.002452 


± 


0.004579 


± 


0.000700 


1.04 - 1.12 


0.251782 


± 


0.113387 


± 


0.068974 


0.278853 


± 


0.108155 


± 


0.064602 






0 






1.12 - 1.20 


0.001304 


± 


0.029359 


± 


0.000143 


0.001304 


± 


0.029359 


± 


0.000143 






0 
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TABLE IX: The diphoton production cross section differential in the cosine of the polar angle in Collins-Soper frame. The first 
error in the cross section is statistical and the second systematic. 



cos 9 bin 


Cross section without cut 


Cross section for Pt<M 


Cross section for P 


T>Af 






[pb] 






[pb] 






[pb] 






-1.20 - -1.12 




0 






0 






0 






-1.12 - -1.04 




0 






0 






0 






-1.04 - -0.96 


0.696711 


±0.146723 


±0.155742 


0.000604 


± 0.012602 ± 0.000066 


0.679356 


± 


0.141129 


± 


0.153197 


-0.96 - -0.88 


0.593058 


±0.160701 


±0.183348 


0.043900 


± 0.032514 ± 0.018749 


0.534910 


± 


0.155592 


± 


0.157684 


-0.88 - -0.80 


0.901409 


±0.165027 


± 0.221503 


0.214033 


± 0.066821 ± 0.046916 


0.679008 


± 


0.152504 


± 


0.179300 


-0.80 - -0.72 


0.964840 


±0.155297 


± 0.234473 


0.686284 


±0.131504 ±0.167953 


0.394443 


± 


0.116872 


± 


0.094604 


-0.72 - -0.64 


2.322500 


± 0.330788 


± 0.674056 


1.895647 


± 0.321399 ± 0.603007 


0.702685 


± 


0.163956 


± 


0.152732 


-0.64 - -0.56 


3.911811 


± 0.353577 


± 0.984306 


3.971193 


± 0.364919 ± 0.976839 


0.308842 


± 


0.140763 


± 


0.123241 


-0.56 - -0.48 


4.225161 


± 0.385762 


± 1.118926 


4.381881 


±0.405094 ± 1.105399 


0.225221 


± 


0.114141 


± 


0.141418 


-0.48 - -0.40 


5.878916 


± 0.500103 


± 1.841742 


5.998004 


± 0.523626 ± 1.862583 


0.328751 


± 


0.123479 


± 


0.118609 


-0.40 - -0.32 


8.447426 


± 0.558218 


± 2.115941 


8.476556 


±0.589264 ±2.168045 


0.611133 


± 


0.126341 


± 


0.129931 


-0.32 - -0.24 


8.006947 


± 0.563629 


± 2.281457 


8.205863 


± 0.595600 ± 2.316283 


0.424552 


± 


0.122576 


± 


0.141255 


-0.24 - -0.16 


9.621518 


± 0.610748 


± 2.601253 


10.100368 ± 0.646808 ± 2.669680 


0.352094 


± 


0.144908 


± 


0.177738 


-0.16 - -0.08 


11.031545 ± 0.647516 ± 3.074353 


11.400319 ± 0.685597 ± 3.182760 


0.544025 


± 


0.146581 


± 


0.146555 


-0.08 - 0 


11.701358 ± 0.652059 ± 3.229385 


12.252276 ± 0.687876 ± 3.313652 


0.409344 


± 


0.154987 


± 


0.199857 


0 - 0.08 


10.744885 ± 0.639661 ± 3.180184 


10.786094 ± 0.672314 ± 3.242347 


0.808290 


± 


0.156327 


± 


0.195042 


0.08 - 0.16 


10.215599 ± 0.620273 ± 2.867342 


10.599528 ± 0.657294 ± 2.971375 


0.475211 


± 


0.141936 


± 


0.137377 


0.16 - 0.24 


8.894859 


± 0.591847 


± 2.681245 


9.066109 


±0.621564 ±2.735215 


0.593223 


± 


0.163517 


± 


0.176269 


0.24 - 0.32 


8.996937 


± 0.573034 


± 2.308823 


9.083449 


± 0.602347 ± 2.317747 


0.814932 


± 


0.184770 


± 


0.227355 


0.32 - 0.40 


8.131862 


± 0.550158 


± 2.082427 


8.160151 


±0.578638 ±2.111536 


0.687574 


± 


0.151101 


± 


0.156505 


0.40 - 0.48 


7.136364 


± 0.535862 


± 1.790615 


7.337457 


± 0.558791 ± 1.797939 


0.363666 


± 


0.150402 


± 


0.157323 


0.48 - 0.56 


6.314780 


± 0.443943 


± 1.479071 


6.414981 


± 0.468845 ± 1.498226 


0.642568 


± 


0.164994 


± 


0.156587 


0.56 - 0.64 


4.206630 


± 0.370591 


± 1.022774 


4.311754 


± 0.395680 ± 1.044270 


0.438339 


± 


0.137911 


± 


0.112795 


0.64 - 0.72 


2.507368 


± 0.287244 


± 0.605107 


2.372216 


± 0.277399 ± 0.544206 


0.242620 


± 


0.118928 


± 


0.127240 


0.72 - 0.80 


0.723595 


±0.163114 


± 0.254653 


0.508954 


±0.128116 ±0.162397 


0.253995 


± 


0.126498 


± 


0.118790 


0.80 - 0.88 


0.553629 


±0.131739 


±0.141735 


0.117661 


± 0.060000 ± 0.032288 


0.422761 


± 


0.113453 


± 


0.106068 


0.88 - 0.96 


0.486210 


±0.125296 


±0.141378 


0.046081 


± 0.036977 ± 0.015266 


0.420675 


± 


0.114459 


± 


0.120587 


0.96 - 1.04 


0.389773 


±0.141480 


±0.142265 


0.014859 


± 0.023566 ± 0.010889 


0.372354 


± 


0.138736 


± 


0.128626 


1.04 - 1.12 




0 






0 






0 






1.12 - 1.20 




0 






0 






0 
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TABLE X: The diphoton production cross section differential in the sub-leading to leading photon transverse momentum ratio. 
The first error in the cross section is statistical and the second systematic. 



2 bin 


Cross section without cut 
[Pb] 


Cross section for Pt<M 
[pb] 


Cross section for Pt>M 
[pb] 


0 - 0.033 


0 


0 


0 


0.033 - 0.067 


0.001845 ± 0.039308 ± 0.000183 


0 


0.001845 ± 0.039308 ± 0.000183 


0.067 - 0.100 


0.064683 ± 0.052631 ± 0.011964 


0 


0.064683 ± 0.052631 ± 0.011541 



0.100 - 0.133 0.040160 ± 0.027228 ± 0.018852 0.014680 ± 0.032433 ± 0.004628 0.059705 ± 0.042528 ± 0.042287 

0.133 - 0.167 0.208681 ± 0.066194 ± 0.042837 0.011290 ± 0.006518 ± 0.002107 0.219178 ± 0.080746 ± 0.052940 

0.167 - 0.200 0.581526 ± 0.144204 ± 0.158584 0.130976 ± 0.080148 ± 0.047058 0.421745 ± 0.112886 ± 0.095122 

0.200 - 0.233 0.687775 ± 0. 179968 ± 0. 138709 0.220290 ± 0.099866 ± 0.041656 0.624676 ± 0. 198432 ± 0. 138087 

0.233 - 0.267 0.882504 ± 0.219745 ± 0.200850 0.311307 ± 0.158139 ± 0.120736 0.662516 ± 0.169350 ± 0.152908 

0.267 - 0.300 1.254394 ± 0.290529 ± 0.272821 0.891524 ± 0.232063 ± 0.173525 0.380905 ± 0.188202 ± 0.125645 

0.300 - 0.333 0.742542 ± 0.342084 ± 0.399154 0.427234 ± 0.266842 ± 0.247371 0.362702 ± 0.241421 ± 0.173907 

0.333 - 0.367 0.633573 ± 0.347426 ± 0.629057 0.465991 ± 0.313708 ± 0.505393 0.171920 ± 0.157777 ± 0.130833 

0.367 - 0.400 2.120037 ± 0.452224 ± 0.852624 1.715373 ± 0.410393 ± 0.695311 0.360338 ± 0.169736 ± 0.140800 

0.400 - 0.433 3.736916 ± 0.590222 ± 1.123369 2.770300 ± 0.524436 ± 0.913390 1.008748 ± 0.280388 ± 0.241464 

0.433 - 0.467 3.631943 ± 0.628078 ± 1.112420 2.889061 ± 0.570550 ± 0.962623 0.784327 ± 0.273687 ± 0.172160 

0.467 - 0.500 4.521544 ± 0.752621 ± 1.788419 3.644498 ± 0.692765 ± 1.516931 0.862816 ± 0.286054 ± 0.268761 

0.500 - 0.533 4.504285 ± 0.765007 ± 2.208000 3.370450 ± 0.701471 ± 2.016244 1.208186 ± 0.316267 ± 0.270012 

0.533 - 0.567 5.323637 ± 0.834206 ± 2.560178 4.458064 ± 0.775756 ± 2.247832 0.852710 ± 0.299524 ± 0.308661 

0.567 - 0.600 7.041574 ± 0.923763 ± 2.652894 6.228433 ± 0.870212 ± 2.309346 0.758477 ± 0.289013 ± 0.321800 

0.600 - 0.633 8.403337 ± 1.002178 ± 3.176046 7.326917 ± 0.951026 ± 2.835386 1.050768 ± 0.308693 ± 0.336094 

0.633 - 0.667 10.720245 ± 1.130024 ± 3.698758 9.382812 ± 1.061512 ± 3.277018 1.561377 ± 0.447592 ± 0.490694 

0.667 - 0.700 12.252908 ± 1.215009 ± 4.376297 11.087163 ± 1.154250 ± 4.032151 1.190425 ± 0.387427 ± 0.352079 

0.700 - 0.733 16.415794 ± 1.325045 ± 5.439192 14.773602 ± 1.265468 ± 4.943280 1.709796 ± 0.404441 ± 0.516840 

0.733 - 0.767 24.203773 ± 1.475829 ± 6.110600 21.437672 ± 1.403616 ± 5.494815 2.989357 ± 0.484161 ± 0.680585 

0.767 - 0.800 25.442360 ± 1.568095 ± 6.926982 23.305996 ± 1.500378 ± 6.332095 2.170811 ± 0.461754 ± 0.605827 

0.800 - 0.833 26.262505 ± 1.639982 ± 7.542237 24.221176 ± 1.570060 ± 6.897027 2.164903 ± 0.503890 ± 0.696416 

0.833 - 0.867 33.239754 ± 1.758692 ± 8.786769 31.182541 ± 1.696572 ± 8.184977 2.019600 ± 0.455356 ± 0.595763 

0.867 - 0.900 43.423592 ± 1.926264 ± 10.230121 41.122257 ± 1.860713 ± 9.587927 2.478693 ± 0.538976 ± 0.732571 

0.900 - 0.933 41.876434 ± 1.940630 ± 9.835004 40.736725 ± 1.885258 ± 9.372191 1.262268 ± 0.497901 ± 0.732053 

0.933 - 0.967 46.414635 ± 2.069529 ± 10.663040 44.152313 ± 2.017412 ± 10.112022 2.116757 ± 0.431365 ± 0.515146 
0.967 - 1.000 43.381981 ± 2.094066 ± 10.489670 41.883457 ± 2.042450 ± 9.953175 1.217864 ± 0.363204 ± 0.483372 





